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Measurements of high frequency temperature and- velocity 
fluctuations were made in the atmospheric boundary layer. 
These data were analyzed using analog correlating techniques 
to obtain the structure function rather than spectral methods 
to obtain the spectral densities. The autocorrelation func- 
tion of temperature fluctuations at 4 meters and 2 meters was 
computed as was the autocorrelation function of downstream 
velocity fluctuations at 4 meters. The cross-correlation 
function of temperature and downstream velocity fluctuations 
at 4 meters was also examined. These results supported the 
Kolmogorov ideas at high wave numbers . There is apparently 
very little correlation between downstream velocity fluctua- 
tions and temperature fluctuations. This indicates the inde- 
pendence of the two sensing systems and that temperature 
fluctuations are related primarily to vertical air movement 
rather than random fluctuations in the horizontal air flow. 
This aspect is further clarified by examining the temperature 
fluctuation records at two heights. 
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I. 



INTRODUCTION 



The need to develop methods of analyzing high frequency 
turbulence signals is apparent. More emphasis is being given 
to the collection of high frequency turbulence data both in 
the atmosphere and in the ocean. This is partly due to the 
fact that technological developments have made such measure- 
ments feasible and, equally important, due to an increasing 
need for such data in directly estimating fundamental turbu- 
lent parameters such as energy and scalar dissipation rates. 

Analysis of high frequency turbulence data presents some 
difficulties not encountered in low frequency analysis. Ex- 
tremely large amounts of information are contained in rela- 
tively short records. In many cases this means the data must 
be analyzed in real time since the compression of the record 
by playback speed-up is not usually possible. If analog data 
are digitized, sampling rates must be very high which limits 
the length of data record that may be transferred to digital 
tape. High frequency turbulence signals tend to be intermit- 
tent so that areas of strong correlation, zero correlation or 
low correlation may be found over relatively short intervals 
of two signals. Furthermore, high frequency signals are small, 
often not much greater than system noise level. 

These are some of the problems investigated in this the- 
sis. First, how does one go about analyzing such data; and 
second, once analyzed how are results to be interpreted? 
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Correlation and spectral analyses have found wide appli- 
cation in engineering and physical investigations for some 
time. More recently the techniques and experiences of these 
established sciences have been carried over into investigating 
geophysical phenomena. Perhaps the outstanding example in 
meteorology and oceanography of the application of these tech- 
niques is in micrometeorology and air-sea interaction. The 
correlation between temperature fluctuations and vertical 
velocity fluctuations leads to a direct estimate of sensible 
heat flux, one of the major terms in thermally -driven models 
of the atmosphere. The correlation between downstream and 
vertical velocity fluctuations leads to a direct estimate of 
the Reynolds stress, or momentum flux, which plays a dominant 
role not only in wind-generated wave theory but also in large- 
scale wind-driven ocean circulation. Spectral analyses, com- 
bined with turbulent theories, have allowed direct estimates 
of the rate of turbulent kinetic energy dissipation and the 
rate of scalar dissipation of such quantities as temperature 
and humidity fluctuations. Such estimates require measure- 
ments of turbulent fluctuations to the highest possible res- 
olution of sensors and electronics. The distribution of these 
quantities in space and time leads to an increase of under- 
standing of turbulent processes. Whereas clearly the goal is 
the understanding of the large scale dynamics of the atmos- 
phere and ocean, such understanding cannot be achieved until 
the very small processes and the non-linear ways in which they 
enter into large processes are understood. 
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One application of studies in this field which is of 
special interest to the military is the possible improvement 
of radio communications based upon an understanding of atmos- 
pheric turbulence. Electromagnetic waves are refracted and 
scattered in the atmosphere. The nature of this scattering 
and refraction is determined in large part by the local re- 
fractive index of air, which in turn is dependent upon both 
the large and small scale temperature fluctuations present. 
Changes in the refractive index can change the path of radio 
waves and seriously disrupt both long and short range radio 
communications. Therefore, increased understanding of tem- 
perature fluctuations and the statistics of temperature fluc- 
tuation distribution, even at very small scales, can lead to 
an improvement of communications by providing the necessary 
information to allow correct action in changing location of 
transmitter or receiver, changing frequency, or changing the 
method of propagation to suit atmospheric conditions. 



8 



II. DESCRIPTION OF DATA 



Based upon the desire for information concerniftg the dis- 
tribution of temperature within a turbulent velocity field 
and the possible correlation in such a field between temper- 
ature and velocity fluctuations, it was decided to measure 
simultaneous temperature and velocity signals with instrumen- 
tation as sensitive as possible. The difficulty of accurate- 
ly measuring temperature fluctuations in a velocity field led 
to the development of a platinum wire resistance thermometer 
by Boston (1970) . This instrument was used to measure tem- 
perature fluctuations and a hot-wire anemometer was used to 
measure velocity fluctuations . Measurements of temperature 
were taken at 4 and 2 meters above the ground, and velocity 
was recorded at 4 meters above the ground. 

A. FIELD SITE 

The data examined in this study were collected at Boundary 
Bay, British Columbia (Figure 1) . This location was the third 
of three sites considered. The first, Spanish Banks, British 
Columbia was not entirely satisfactory because temperature 
signals tended to be small in the stable air over the ocean. 
The second, Hanford, Washington, was not satisfactory because 
of the varied meteorological conditions encountered in air 
over a desert. The third. Boundary Bay, was a compromise be- 
tween the first two. The air was slightly unstable as it 
blew across the beach off the ocean. The site is a tidal 
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Figure 1. Location of Field 
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mud-flat, which is completely covered at high tide and is 
bare for about 2.5 km seaward at low tide. The site has been 
described in detail by Dobson (1969) . 

B. EXPERIMENTAL EQUIPMENT AND SENSORS 

All instruments were supported from a portable mast at 
heights of either 4 or 2 meters. The signals from the sensors 
were carried by cables to a panel truck which served as the 
housing for the recording equipment. The truck was 30 meters 
downwind from the mast. 

1. Platinum Resistance Thermometer 

The sensor of temperature fluctuations was a platinum 
resistance thermometer consisting of a wire 0.25y in diameter 
and 0.25 mm in length. The wire formed one leg of a Wheat- 
stone bridge. A square wave of constant amplitude is applied 
to the bridge via an 80 kHz multivibrator and switching tran- 
sistor. A synchronous detector is used on the bridge output. 
D.C. amplification is performed by an operational amplifier 
with fixed incremental voltage gains. Its frequency response 
is flat from D.C. to approximately 10 kHz. The platinum re- 
sistance thermometer electronics were placed at the foot of 
the mast . 

2 . Hot-wire Anemometer 

A constant temperature hot-wire anemometer system 
(DISA model 55D05) was used. The unit was battery operated 
and used a 1:1 bridge ratio made with a compensating cable. 
The hot-wire probe consisted of a platinum wire 5y in diam- 
eter and 1.0 mm in length. The hot-wire was mounted 
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approximately 5 cm above the temperature sensor. Elec- 
tronically this system can handle signals in the frequency 
range from D.C. to 100 kHz. Both the time response and wave 
number limitation of the wire allow signals up to 1 kHz to 
be measured without attenuation at wind speeds of 4 to 
5 m/sec. 



3. Cup Anemometer and Quartz Thermometers 

In addition to the above instruments, a cup anemom- 
eter and two quartz thermometers were used in the experimen- 
tal setup. The cup anemometer, mounted on the mast at 4 
meters, was used to measure mean wind speed in order to pro- 
vide (1) an in-field check on the hot-wire anemometer, and 
(2) a mean wind speed, U (cm/sec), to apply in Taylor's 
hypothesis in the form 



k 



2 TT f 
U 



( 1 ) 



where k is the radian wave number 



(cm“^) 



and f is frequency 



(Hz) . 



The quartz thermometers were used to get an indica- 
tion of the relative stability of the atmosphere. This data 
would be important over a 2 4 -hour period, but not over a 
period as short as 8 minutes, which was the length of the 
recorded data used in this thesis. 

4 . Recording System 

The signals from the platinum resistance thermometer 
and hot-wire and cup anemometers were recorded on separate FM 
(frequency modulated) channels of a magnetic tape recorder 
(Sangamo model 3562) run at 60 inches per second. The high 
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frequency cut-off of the tape recorder was 20 kHz . The 
quartz thermometer digital meter reading was read onto the 
voice edge track of the tape. The recording system is out- 
lined in Figure 2 . 

Since both the noise of the system and the maximum 
frequency of interest of the temperature signal were uncer- 
tain at the time of the experiment, recordings were made both 
without the low pass filter (Krohn-Hite model 3340) and with 
it at various settings to attenuate very high frequency noise. 
Data chosen to be analyzed in this thesis were recorded with a 
low pass filter cut-off at 10 kHz (-3 db down frequency) . 

The signals from the platinum resistance thermometer 
(PRT) and the hot-wire anemometer were treated similarly. 

Each was recorded in two different ways chosen to optimize 
the incompatible requirements of noise performance and fre- 
quency response. One channel was unfiltered (L.F. GAIN) and 
the other was differentiated. The purpose of the differenti- 
ation circuit was to improve the signal-to-noise ratio at high 
frequencies so that the very small scale fluctuations could 
be more readily examined. L.F. GAIN provided gain not only 
for the direct signal but also for the signal prior to dif- 
ferentiation. The high frequency gain control (H.F. GAIN) 
determined the frequency at which there was unity gain in the 
differentiating circuit. H.F. GAIN was adjusted to provide 
optimvim gain (prewhitening) of the frequencies of interest. 

The differentiating circuits were matched linear 
phase shift devices made from Philbrick P85AU operational 
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Figure 2. Block Diagram of Recording System 
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amplifiers. The circuit had a 6 db/octave gain up to 10 kHz 
(depending on setting) and then fell at 6 db/octave. 

Recording at 60 in/s with FM electronics ensured the 
best signal-to-noise ratio and provided sufficient’ frequency 
response (up to 20 kHz) to record faithfully the sharp rises 
and falls of the differentiated signals. It has the disad- 
vantage of limiting the maximum length of a given record to 
8 to 12 minutes, depending on the length of the tape. 

Specifically in the case of. the data accumulated for 
analysis in this thesis, measurements of temperature and 
velocity fluctuations took place at Boundary Bay on the after- 
noon of August 8, 1969. Measurements were recorded on mag- 
netic tape No. 20 3(1) . Record length was 8 minutes commencing 
at 1427 PDT. A platinum resistance thermometer was located 
on the mast 4 meters above the ground. The temperature signal 
was passed through the PRT electronics where the signal was 
amplified, then through the Krohn-Hite low pass maximum flat 
filter with the cut-off frequency set at 10 kHz. The output 
of the Krohn-Hite filter was fed to the Philbrick amplifier 
manifold where both the low frequency and differentiated sig- 
nals were amplified again. The low frequency, undifferen- 
tiated signal was then recorded on tape channel 5. The high 
frequency differentiated signal was recorded on channel 7. 

The DISA hot-wire anemometer also was located on the 
mast 4 meters above the ground. The velocity fluctuation 
signal from this sensor was fed through a Philbrick amplifier 
manifold for signal gain. The low frequency, undifferentiated 



15 



signal was then recorded on tape channel 2. The high fre- 
quency, differentiated signal was recorded on channel 4. 

In addition, at 2 meters elevation there was another 
platinum resistance thermometer. The temperature fluctuation 
signal from this sensor was fed through and amplified by 
separate PRT electronics, and then passed to a Krohn-Hite low 
pass maximum flat filter with a cut-off frequency setting of 
10 kHz. The output of the Krohn-Hite filter was then again 
amplified by the Philbrick amplifiers. The low frequency, 
undifferentiated signal was then recorded on channel 1, and 
the high frequency, differentiated signal on channel 3. 

In the case of all three velocity and temperature 
measurements, each signal was amplified by the PRT electronics 
and the Philbrick operational amplifiers to such a degree that 
complete advantage was taken of the full dynamic range of the 
tape recorder. 
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III. THEORETICAL RELATIONS 



Correlation analysis is a method of measuring the inter- 
dependence of variables associated with various processes, 
often of a random or periodic nature. Simultaneous time 
functions, such as turbulence data measured at a fixed point 
in space, can be described by correlation functions. These 
functions may be of the dependence of a variable at one time 
with itself at another time (called the autocorrelation func- 
tion) or the dependence of a variable with another variable 
(called the cross correlation function) . 

A. AUTOCORRELATION AND CROSS CORRELATION FUNCTIONS 

The autocorrelation function of a time series g^(t) at 
the time t and (t + x) will be indicated by . It is 

obtained by averaging the products of two values separated by 
a lag time, x, over an observation time, T. The resulting 
average product is a function of x: strictly 



At X = 0 , R. . (x) has its maximum value which is related to 

t XI ' 

the variance of the signal by 



T 




o 



( 2 ) 



= g^(t) g^(t+x) 



( 3 ) 




( 4 ) 
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for all T. Clearly 



Rii(T^) = • (5) 

The form of Rj_jL (^ ) will be characteristic of the original 
signal g^^ (t) . 

The autocorrelation function is related to the energy 
spectrum (sometimes called the power spectrum) as 

follows ; 



CO 

R^j^(x) = J' E^^(f) • cos 2irfT • df • (6) 

o 

and the inverse is: 

oo 

E^jl^(f) = 4 J" * cos 2iTfT * (7) 

o 

which merely states that R. . (t) and E. . (f) are Fourier cosine 

j 11 11 

transforms of each other. One can be found from the other by 
integration. 

The theoretical relations involved in performing cross 
correlation of a time series are the same as those for the 
autocorrelation function described here with the exception 
that two different signals are being correlated rather than 
the same signal. That is, instead of correlating gj^(t) with 
itself, gji^ (t) is correlated with gj (t) . 



18 



B. KOLMOGOROV'S POWER LAWS AND THE STRUCTURE FUNCTION 

Kolmogorov (1941) postulated that high frequency velocity 
fluctuations should be isotropic even though they are gener- 
ated from and embedded in anisotropic eddies of low frequen- 
cies. If there is a significant region of the spectrum at 
high frequencies in which there is no energy feeding, then 
since there is an energy cascade from lower frequencies by 
inertial transfer, with the pressure forces tending to trans- 
fer energy among components, it would be reasonable to suppose 
that these smaller scales would not exhibit the anisotropy 
of the lov7 wave numbers and would tend toward isotropy. This 
is known as local isotropy, and has been supported by experi- 
mental evidence of Batchelor (1953). 

If the Reynolds numbers are high enough, turbulence will 
adjust through inertial transfer and viscous dissipation 
until a statistical steady-state, independent of the initial 
conditions, is obtained. This steady-state consists of energy 
transfer into the high frequency range that is equal to energy 
lost through viscous dissipation. 

Kolmogorov hypothesized that at high Reynolds numbers in 
any turbulence, the small-scale, high frequency eddies are 
isotropic and in universal statistical steady-state. Their 
properties should depend only on the rate of dissipation of 
energy (e) and on the kinematic viscosity (v) , and not direct- 
ly on time. Therefore motion associated with the steady- 
state range is uniquely determined by the parameters e and 
V . Based on dimensional analysis, Kolmogorov's conclusion 
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about the form of the three-dimensional energy spectrum v/as 



E.i(k) 



T, /k . 

= F(^) 

s 



(8) 



(ev5)V4 

where F is a universal function, k is the wave number, is 
the Kolmogorov wave number, and l/k^ is the Kolmogorov 
microscale . 

This has been verified quite conclusively by Grant, 
Stewart, and Moilliet (1962) and Pond, Stewart, and Burling 
(1963) . From their completely separate measurements in’ a 
tidal channel and in the atmosphere, all results agreed with 
the above functional relation over 4.5 decades of wave 
nvunbers . 

Kolmogorov's second hypothesis concerned the inertial sub- 
range. If the Reynolds number is large enough, say, in the 
order of 10^, the energy-containing eddies and the dissipa- 
tion wave number regions should be well separated in wave 
number space. There should, therefore, exist a range of wave 
numbers on the low wave number side of the dissipation range 
but at higher wave numbers than the energy-containing wave 
numbers in which universal statistical steady-state holds, 
but at which there is negligible dissipation of energy. The 
dominant feature would be an energy cascade by inertial forces 
through this inertial subrange. Since viscosity is not im- 
portant in this wave number range, statistical properties are 
determined by the rate of dissipation only. Thus, in the 
inertial subrange by Kolmogorov's theory, (k) should have 

the following functional relation: 
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(9) 



XI 



(k) 



Ke2/3k-V3 



where K is the Kolmogorov constant. This equation is empir- 
ical and is obtained by dimensional analysis. 

The results presented in this thesis represent structure 
functions of the temperature and velocity fluctuations in the 
atmosphere. The structure function, D(r), depends on veloc- 
ity or temperature differences between two points in the 
turbulent flow, i.e., for velocity; 

D.. (r) = TTT2 (10) 

XL [u(x+r) - u(x)] 



where x is a point in the turbulent field and r is the 
separation distance. 

If the turbulence is homogeneous, the structure function 
is relatea to the correlation function, R^^(r), by 

D^^(r) = 2[R^.(o) - R..(r)]. (11) 

The inverse transform relation between the correlation 
fxmction and the energy spectrum, (k) , is 



R,i(r) = 



00 

JfJ 



Eii (k) • COS kr • dk. (12) 



The relation between the structure function and the spectrum 
is therefore 



Di.(r) 



00 

-■/// 



(1 - cos kr) • 



E. . (k) 

XL 



dk 



( 13 ) 



For conditions of isotropy, equation (13) reduces to 
(Strohbehn 1968); 
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dk . 



(14) 



D..(r) = Sir/* (1 - . E..(k) . k2 . 

XI J kr XX 

o 

It can be seen from equation (14) that the structure func- 
tion, as the spectrum (k) , should also take on a charac- 
teristic power law behavior in isotropic turbulence. From 
dimensional reasoning and the assumption of dependence on 
e and r, this power law relation is 

Di^(r) = C r^/^ , (15) 

where C is a constant of proportionality. 

Similar arguments may be applied to temperature fluc- 
tuations to yield expressions for the temperature spectrum and 
temperature structure function with the same power law be- 
havior. For isotropic turbulence, the temperature structure 
function can be related to the temperature autocorrelation 
function by 



(Tf - ^ 2T^(1 - R^^) (16) 

where R. . is the autocorrelation function. Thus 

XX 

(1 - R^^) a r^/^ . (17) 

This relation is expressed in terms of the separation 
distance, r, whereas the correlation analysis was made by 
varying lag time, t. Assuming homogeneity of the turbulence 
field, which seems entirely justified for the observations 
made in this study, the correlation in space should be equiv- 
alent to the correlation in time. The relation between r and 
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T is 



T 



r 

U 



(18) 



where T is the lag time and U is the mean wind speed. This 
is a Variation of Taylor's frozen turbulence hypothesis, 
which .is most often used to transform frequency (f) (time) 
spectra to wave number (k) (space) spectra; 

k = . (19) 

From (18) and (19) , one may relate t to frequency since 

( 20 ) 
( 21 ) 



1 k . 
r 



So f = 



2 7T T • 
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IV. ANALYSIS PROCEDURES 



Because of the need to develop methods of analyzing high 
frequency turbulence signals, it was decided to use two dis- 
tinct analysis procedures. An attempt was first made to 
perform spectral analysis on a 6-minute record of data by 
digital procedures. This was accomplished by means of analog- 
to-digital conversion (XDS 9300/CI 5000 Hybrid Computer System) 
followed by IBM-360/67 digital computer spectral analysis. 

This procedure met with only limited success and is explained 
in Appendix A. The second procedure consisted of a correla- 
tion analysis of the analog temperature and velocity fluctua- 
tion signals by a DISA Analog Correlator. The results of this 
second method of analysis comprise the major results of this 
thesis . 

A. SELECTION OF DATA 

Since the objectives of this thesis were concerned with 
the examination of high frequency turbulence signals, it was 
not necessary to examine long sections of data. Thus it was 
decided that a record of about one minute length would be 
sufficient. Actually the exact record length turned out to 
be 50 seconds because of the scale used in recording the out- 
put of the DISA Analog Correlator. 

The primary criterion for selecting a section of data for 
analysis was stationarity of velocity signals. In addition 
it was desirable that the selected section contain 
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temperature fluctuations at both 4 and 2 meters elevation that 
could be considered typical of the 8-minute record. There- 
fore, with these ideas in mind, the data record selected was 
a 50-second section of tape 203(1) commencing 22 seconds 
after the beginning of the measurements. This section ap- 
pears in Figure 3. Note that high frequency turbulence inten- 
sity of both the direct and differentiated velocity signals 
remains reasonably uniform over the length of this section. 
Figure 4 shows the direct signals only. The direct velocity 
signal shows that some low frequency oscillations are present, 
but it is not difficult to see that the velocity signal devi- 
ates very little from the mean wind speed of 3.0 m/sec. 

A variety of temperature fluctuations were encountered 
throughout the 8-minute record. Some sections showed a signal 
that tended to be one-sided but with the small scale fluctua- 
tions remaining uniformly active throughout. The one-sided 
feature is a shape which is often characteristic of tempera- 
ture signals; a gradually rising leading edge and a sharp 
drop-off of the trailing edge, which gives the record a saw- 
tooth pattern. In general the signal exhibited intermittent 
quiet and active periods. The 50-second record selected at- 
tempted to incorporate as many of these fluctuation character- 
istics as possible. These characteristics will be interpreted 
in more detail in view of the scope and objectives of this 
thesis in Section V. 
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Figure 3. Data Record Selected for Analysis 

(Direct and Differentiated Signals) 
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Figure 4. Data Record Selected for Analysis (Direct Signals) 
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B. CORRELATION EQUIPMENT AND PROCEDURES 

The correlation analysis was performed using a DISA Type 
55D70 Analog Correlator System. In support of the DISA cor- 
relator were a DISA Type 55D75 Time Delay Unit and a DISA 
Type 52B01 Sweep Drive Unit. The sweep drive unit was used 
to provide a fractional breakdown of the delay range setting 
of the time delay unit. This feature permitted the use of 
convenient fractional lag times within the range of the set- 
ting on the time delay unit. In addition the sweep drive 
unit provided control for the Hewlett/Packard Moseley 135 X-Y 
Recorder which was used to record the correlator output. 

The signals being correlated were reproduced on the Sangamo 
3562 tape recorder and fed to the time delay unit. The speed 
of the tape recorder for all correlation measurements was 
60 in/s (real time) . The delay range setting of the time 
delay unit was 0-10 ms when correlating undifferentiated 
signals and 0-1 ms when autocorrelating differentiated signals. 

The signal was then fed through the DISA correlator and 
sweep drive unit. The integrating time constant setting used 
in all correlations was 10 seconds. This provided enough 
signal smoothing to recognize the large fluctuations in cor- 
relation, while also smoothing the signal of small scale 
fluctuations irrelevant to the fixed-lag correlation study. 

The lag times set into the sweep drive unit varied with the 
correlation of different signals. The sweep drive unit had 
the capability of reducing the delay range setting into frac- 
tions of 1/500 . Thus different settings of lag time could 
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be chosen to adequately cover the correlation range from 
-1.0 to 1.0. Some specific settings for selected autocor- 
relations are given in Section V. 
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V. RESULTS 



The results of this study are significant in light of its 
objectives, which were to (1) investigate methods of examin- 
ing universal properties of high frequency turbulence signals 
and (2) to specifically examine the relation between tempera- 
ture fluctuations and downstream velocity fluctuations in the 
atmospheric boundary layer. A visual analysis of the tempera- 
ture fluctuations, especially at 4 meters elevation, provided 
substantial support of the work of Taylor (1958) and Priestly 
(1959) . The correlation analysis of both direct and differen- 
tiated signals at 2 and 4 meters elevation obeyed correlation 
analysis theory and also agreed with Kolmogorov's -5/3 power 
law relation for the inertial subrange when the data were 
analyzed by the structure function. Finally, a lack of cor- 
relation between velocity and temperature fluctuations at the 
same height indicated the independence of the two sensing sys- 
tems and therefore the fact that their design is good. 

A. EVIDENCE OF ORGANIZED MOTION IN THE FLOW 

A consistent characteristic of the temperature data was 
the occurrence of "pulses" of large amplitude. These pulses 
occurred nearly simultaneously at both 2 and 4 meters eleva- 
tion. Their occurrence was intermittent and always repre- 
sented an increase in temperature. They occurred within 
quiescent periods during which the recorded trace appeared to 
have a preferred mean value. 
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A selected portion of a strip chart record exhibiting 
these pulses appears in Figure 5. The pulses (e.g., point A 
in Figure 5) have a triangular shape which is distorted to 
give a saw-tooth pattern. Taylor (1958), Webb (1965) and 
Priestly (1959) have reported similar occurrences in data of 
these kind. 

Taylor (1958) argues convincingly that the pulses are due 
to organized thermal structures, "plumes," embedded within 
the turbulent flow. Taylor reached his conclusion by perform- 
ing correlation analyses between temperature fluctuations at 
several levels. The time lag for maximum correlations was in 
agreement with an expected tilt of these plumes due to the 
shear flov;. A schematic of the plume in a shear flow appears 
in Figure 6 . 

A time lag between the occurrence of the pulse at the 2 
and 4 meter level is evident at Point A in Figure 5. This 
lag indicates that the plume was swept past the 4 meter level 
before the 2 meter level. Because the sensors were mounted in 
a plum-line arrangement, these results are in agreement with 
the concept that plumes, under the influence of the wind 
shear, are responsible for the observed pulses. 

The saw-tooth characteristic of the pulses is, perhaps, 
due to the plumes having a well-defined trailing edge, but a 
leading edge that is altered by mixing or entrainment. Mixing 
would lower the temperature in the leading portion of the 
plume giving the saw-tooth pattern . 
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Figure 5. Strip Chart Record of Temperature Pulses Exhibiting Plume Properties 
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Figure 6 . Schematic of Thermal Plumes in Shear Flow 



Considering the consistent occurrence of pulses in these 
temperature data in view of Taylor's (1958) more extensive 
study, it is evident that organized motion exists within the 
turbulent flow. By selecting a coordinate system consisting 
only of a trace of the temperature versus time, it is pos- 
sible to obtain a description of the plumes in the complex 
regime. The presence of organized motion does not appear to 
influence the universal properties of the fluctuations at the 
small scales. Properties at high frequencies are examined by 
another coordinate system, the structure function versus time 
lag . 

B. CORRELATION ANALYSIS 

The results of the autocorrelation analysis were consist- 
ent and as expected. The correlation function, Rj^j^ ^ of both 
the direct and differentiated signals decreased as the lag 
time, T, increased. This was true throughout the 50-second 
record. Figure 7 is a typical example of such behavior. 

Here direct temperature fluctuations at 4 meters elevation 
are examined. Note that with t set at 0 ms, the autocorrela- 
tion function exhibits a high value (=:.9 3) and that as the 
lag time is increased, R^^ decreases in value until, with a 
lag time of 10 ms, a negative correlation of mean value -.28 
exists. The fact that a correlation of 1.0 does not result 
from a x setting of 0 ms is attributed to correlation equip- 
ment calibration. When the time delay unit was set for 0 ms, 
in fact a slight time lag was present. 
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Figure 7. Autocorrelation of Direct Temperature Signal 
at 4 Meters for Various Lag Times (Tau) 
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Time (sec) 

Figure 8. Autocorrelation of Differentiated Temperature Signal 
at 4 Meters for Various Lag Times- (Tau) 
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Figure 8 exhibits an example of the autocorrelation 
results for a differentiated signal. In this case, the dif- 
ferentiated temperature at 4 meters elevation was correlated 
with itself. Since in differentiated signals the very high 
frequencies (i.e., very small scale fluctuations) are accen- 
tuated, the lag times used for differentiated signals are one 
order of magnitude smaller than those used in correlation 
analysis of the non-dif ferentiated, direct signals. However, 
the same autocorrelation characteristics appear. When t is 
set at 0, is high, and as the lag time is increased, the 

value of the correlation function decreases. 

In Figures 7 and 8, a definite pattern of autocorrelation 
function fluctuation exists across the entire 50 -second record. 
Times of interesting activity were selected and the structure 
function plotted. 

Figure 9 shows the autocorrelation function of temperature 
fluctuations at 4 meters elevation matched with the 50-second 
analog trace of these same temperature fluctuations. It can 
be seen that the area of the temperature plume at the 36.5 
second point coincides with a definite increase in correlation 
at the same time in the autocorrelation record. 

Figure 10 shows a plot of the cross correlation of direct 
signals of temperature and velocity fluctuations at 4 meters 
elevation. The continuous low values of the correlation 
function, regardless of the lag time, x, indicates that 

these two signals are poorly correlated. This lack of cor- 
relation between simultaneous temperature and velocity 
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Figure 10. Cross Correlation of Direct Temperature and Velocity Signals 
at 4 Meters for Various Lag Times (Tau) 
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fluctuations at the same height indicates that the design of 
the two sensors, the platinum resistance thermometer and the 
hot-wire anemometer, is favorable for measuring such signals. 
A lack of correlation indicates that velocity fluctuations 
are not having a significant effect on the platinum resis- 
tance thermometer; and vice versa, that the temperature fluc- 
tuations do not affect the velocity measurements of the hot- 
wire anemometer. 

C. STRUCTURE FUNCTION 

Figures 11 through 16 show logarithmic plots of the auto- 
correlation structure function, versus lag time, t, 

for all six recorded signals at various times of interest in 
the 50 -second record. In the cases of direct temperature 
fluctuations at 2 and 4 meters elevation and direct velocity 
fluctuations at 4 meters elevation, the slope of the plots 
approximates the +2/3 value predicted by applying 
Kolmogorov's inertial subrange hypothesis to the structure 
function. The slopes of the differentiated signals decrease 
at large values of x. 
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Figure 11. Structure Function Plot of Direct Temperature Signal 
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Figure 12. Structure Function Plot of Direct Temperature Signal 
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T (ms) 

Figure 13. Structure Function Plot of Direct Velocity Signal 
at 4 Meters 
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Figure 14. Structure Function Plot of Differentiated Temperature 
Signal at 2 Meters 
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Figure 15. Structure Function Plot of Differentiated Temperature 
Signal at 4 Meters 
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Figure 15. Structure Function Plot of Differentiated Temperature 
Signal at 4 Meters 
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Figure 16. Structure Function Plot of Differentiated Velocity Signal 
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VI. CONCLUSIONS 



Analysis of the strip chart records of temperature fluc- 
tuations indicated that the data support Taylor's and 
Priestly's description of thermal pliames within a turbulent 
shear flow; and, as such, this study added information on 
these phenomena. 

The data were consistent with predictions by similarity 
theory since both velocity and temperature fluctuations 
assumed a +2/3 power law behavior for the structure function 
in accordance with Kolmogorov's inertial subrange hypotheses. 

The results of this study provide further support that 
the +2/3 power law behavior is valid for temperature fluc- 
tuations (Krechmer, 1952) . Results of previous work in the 
atmospheric boundary layer have been primarily concerned with 
the structure function behavior of velocity fluctuations. At 
first glance, it is surprising that temperature as well as 
velocity fluctuations should exhibit consistent +2/3 slopes 
since the temperature traces showed the presence of intermit- 
tent active and quiet periods associated with thermal plumes. 
These plumes should be associated with anisotropic conditions 
since they imply vertically upward directional preference due 
to buoyancy. Thus, the isotropy deduced from the +2/3 behav- 
ior of the structure function and the anisotropy of the 
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characteristic plumes appear to be a paradox.^ This apparent 
inconsistency can be explained by the fact that in this study 
only the very high frequency temperature fluctuations, which 
are independent of the features of the large scale flow, 
were examined. In order to see the difference in spectra or 
structure functions, the lower frequencies must be examined. 
The examination of these lower frequencies is one obvious 
future extension of this work. 

The favorable design of the two sensors as indicated by 
the lack of correlation between velocity and temperature 
fluctuations at 4 meters elevation suggests that the two sens- 
ing systems used to measure the turbulence in this study are 
compatible and could continue to be used to make further 
measurements . 



A +2/3 slope does not guarantee isotropy; isotropy is 
a sufficient condition for a +2/3 slope, but not a necessary 
condition. For examples, see Weiler and Burling (1967) . 
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APPENDIX A 



DIGITAL SPECTRAL ANALYSIS PROCEDURES 



A significant portion of the research in this study was 
an attempt to digitize the analog signals of temperature and 
velocity fluctuations, and then perform digital spectral 
analysis upon the data by means of a high speed digital com- 
puter. Despite the fact that various problems were encoun- 
tered and only limited success resulted, it is worthwhile to 
discuss the procedures involved in this analysis since they 
could well be applicable to future studies of the relation 
between velocity and temperature fluctuations. 

A. ANALOG TO DIGITAL CONVERSION 

The analog to digital conversion was accomplished by 
means of a XDS 9300/CI 5000 Hybrid Computer system. The XDS 
9300 is a digital computer manufactured by Xerox, Inc. with 
a 32,000 byte storage capacity, and the Cl 5000 is a COMCOR 
analog computer with a 50-amplifier size rating. The temper- 
ature/velocity signal output of the Sangamo 3562 tape recorder 
was connected to a Krohn-Hite model 3340 low pass filter. The 
cut-off frequency settings were 2 kHz for the undifferentiated 
signals and 4 kHz for the differentiated signals. The signals 
were then fed to the analog computer patch panel. The signal 
was amplified in this patch panel to take advantage of the 
full dynamic range of the hybrid computer system. In 
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subsequent troubleshooting, the patch panel was found to be a 
source of substantial noise. It is possible that this was 
one of the problems that precluded good results from this 
procedure . 

From this patch panel, the signals were fed to the XDS 
9300 digital computer. The digitized output of the XDS 9300 
computer was then recorded on 7 -track magnetic tape. Two 
channels of the analog tape recorder were digitized at one 
time. The procedures resulted in digital tape containing 
blocks of data 2048 words in size. Every other word in each 
block was a data sample from one channel. Thus each block of 
data contained 1024 words (data samples) of each channel. 

The sampling rate used for digitizing the undifferentiated 
signal was 2000 samples/second. The differentiated signal 
was sampled at a rate of 4000 samples/second. In the case of 
all signals, a 6-minute section of the 8-minute record was 
digitized. The procedure of digitizing two channels at one 
time was another possible source of error in this process. 

In future work, it is recommended that initially only one 
channel be digitized at a time in order to simplify trouble- 
shooting in later stages of this procedure. Of course, when 
computing cross-spectra, it will be necessary to digitize 
two channels simultaneously. 

For a more specific, detailed description of the general 
analog to digital conversion procedure used in this study, 
including programs, equipment settings, and further recom- 
mendations, see Jones (1971) 
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B. DIGITAL DATA PROCESSING 

The digital data processing was performed on an IBM 
360/67 computer. This system contains a core space of 
762,000 bytes. Supporting equipment includes two compilers, 
two CALCOMP plotters, two printers, and twenty remote units. 
Digital tapes created by the hybrid computer system were 
written in 7-track form (octal representation) . However the 
IBM 360/67 System normally operates in 9 -track form (hexa- 
decimal representation) . Therefore all data derived from the 
analog to digital conversion had to be converted to a 9 -track 
format prior to analysis by the IBM 360/67. Standard sub- 
routines existed for this conversion, but required some modi- 
fications as explained in Jones (1971) . The output of this 
conversion was the digitized data on a 9 -track tape. This 
process was another possible source of error. In trouble- 
shooting this phase of the procedure, a program was written 
to plot the data of the 9 -track tape. Numerous spikes not 
present in the raw analog data were found to be present in 
the 9 -track data. 

The next step in this procedure was the computation of 
the Fourier coefficients for the data. The program FTOR, 
described by Wilson, Boston, and Denner (1969) , was used to 
accomplish this step. The output of the Fourier coefficients 
was recorded on another 9 -track tape. 

The final step in the digital data processing was the 
computation and plotting of the spectra. Two plotting pro- 
grams, SCOR and FCPLOT, described by Wilson et. al . (1969), 
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read the output tape of the FTOR program and plotted various 
spectra. The SCOR program provided a plotted output of 
power spectra and cross spectra, and the FCPLOT program pro- 
duced printer and ink drawings of amplitudes of Fourier coef- 
ficients. Both programs processed a single reel of digital 
tape (1368 blocks of 2048 words each) in less than 10 minutes 
computing time. In this study the SCOR program output of 
spectra and cross spectra for all measured temperature and 
velocity fluctuations exhibited little relation to known power 
spectra behavior for such signals . The FTOR, SCOR, and FCPLOT 
programs have been used successfully in the past by Wilson, 
et . al ♦ (1969) in analyzing such data. This leads to the be- 

lief that the problems encountered and general lack of results 
in this study were probably generated in either the analog to 
digital conversion or the 7-track to 9 -track tape conversion 
processes . 

Detailed information on the FTOR, SCOR, and FCPLOT pro- 
grams can be found in Wilson, et. al. (1969) . Specific de- 
tails of the application of these programs to temperature and 
velocity fluctuations in the atmospheric boundary layer are 
described at length by Jones (1971) . 
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